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YY ABSTRACT 


This report presents the results of a six-degree-of-freedom analytical 
study to examine the stall/near atall/spin characteristics of the slatted 
F-4E aircraft. The basic F-4E and the slatted F-KE are analyzed and compared. 
Primary emphasis was directed toward the stali/near stall characteristics. 
Some analyses of spin characteristics and recovery techniques were nlso 


accomplished. ; 


aes ee 


‘ Departure and post stall characteristics are analyzed in detail. The 
primary points of interest are:\1) types of entry into the stall, (2) the 
effect of lateral-directional control inputs on the depsrture cheracteristics, 

(3) recovery techniques, andf.}) evaluation of five methods that predict the 
point at which the aircraft becomes lateral-directionally «unstable ,) Included 

_-is-a-—new- technique created by McDonneli~Aircraft-Company;> The previous F-4 
spin evaluation developed a dynamic model whereby rotary balance and forced 
oscillation data are utilized to compute the spin and spin entry motions. 
Further development of this dynamic model has been accomplished so that it 
can be used in both the uwistalled and stalled flight regimes. Transition 
between the two factanee ean be smoothly accomplished, ' Based on the investiga- 
tion, recommendations are made relative to the techniques for recovery from a 
departure and for preventing departures through the use of control surfaces, — 
A comparison of the slatted F-lB flight test stall/near stall characteristics 
“with the analytical characteristics are qualitatively evaluated. 

The analytical study of the spin entry and developed spins were performed 

primarily to note any differences in spin characteristics between the basic 
F-lke and the slatted F-NE. The effect of three c.g. positions are studied 


along with the effect of using recovery coatrols in a spin. 
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total angular velocity rad/sec 


component of total aircraft angular velocity parallel 
to velocity vector, positive clockwise looking along 


velocity vector red/sec 
engine rotation rate positive clockwlse looking 

forward rpm, rad/sec 
air density slug/ft3 


MCOONANELL DOUGLAS CORPORATION 
xix 


MCDONNELL AIRCRAFT COMPANY Mpc A2181 


1. Introduction 

In recent years the F-4 has been operated at the extreme of the ma- 
neuvering envelope which requires operating the airplane near stall. As 
@ result an increase in loss of ey and departure/spin incidents have 
occurred. The installation of the two-position leading edge slats, which. 
significantly improves tne maneuvering characteristics of the F-k due to 
improved lateral-directional stability at high angles of attack, should 


reduce the number of such incidents. 


As a contractual requirement, in revrofiting all existing F-4e aircraft 
with loading edge slats, an analytical study was to be made of the stall/ 
near stall/spin characteristics of the pel airplane, The main objectives 
of this study were: 

2) To define the lateral-directional handling 
characteristics of the slatted F-lE. 
2) To make a recommendation on the need for a 
stall/near stall/spin flight test program. 
fo perform this analysis the McDonnell-Douglas Six-Degree-of-Freedom (SDF) 
digital computer program was utilized. Tha reference angle of attack and 
angle of sideslip definitions for this study are presented in Figure 1.1, 
The Euler angle sequence and signconvention used is presented in Figure 1.2. 


Likewise, the thrust vector orientation is presented in Figure 1.3, 


As previously mentioned the addition of leading sdge slats alters ‘the F-HE 
aerodynamic characteristics at high angles of attack. ‘Therefore, a complete 


high angle of avtack wind tunnel test program was undertaken to obtain both 
static and dynamic wind tunnel data. 
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2. Wind Tunnel and Sree Flight Model Test Results For The Basic and 
Slatted F-4E 


2.1 Wind Tunnel Static and Dynamic Investigation - At the start of the 
program studies, a wind tunne). investigation of the F-4 with leading edge 


slats was conducted using various models, 


a) 5% Scale Force & Moment Model 

b) 6.67% Scale Force & Moment Model 

ec) 9.1% Seale Rotary Balance Model 

a) 13% Scale Forced Oscillation Model 
Some additional new data was obtained on the basic F-k in the above tests, 
however, most of the basic F~ wind tunnel data was obtained during the 


previous F-4 analytical spin evaluation study effort (Reference 1). 


2.1.1 Static Force & Moment Wind Tunnel Results - The slatted F~hE 
wind tunnel data obtained were used to generate the static force and ‘moment 
curves used in the analytical study for the slatted F-hE and can be found 
in Appendix.A. A similar set of curves for the basic F-4E can be found 
in Appendix B. These data, which were obtained from Reference 2,' were 
modified where appropriate to provide a better representation of the air- 


plane. 


2.1.1.1 Longitudinal Characteristics ~- The wind tunnel 
pitching moment and normal force curves for the slatted and non~slatted 
configurations are approx-mately the same above 25° wing angle-of-attack 
with some improvement in norma] force due to slats below 20° wing angle- 


of~attack (see Figures 2.1 and 2,2). ‘he axial force above 10° angle-of- 
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attack is more positive due to a suction effect caused by the leading edge 


slats (see Figure 2.3). 


Btabiletor effectiveness data are obtained for 47° to 0° and 0° to 
-21° atabilator settings. Stabilator effectiveness is not affected by the 
addition of leading edge slats. The +79 stabilator setting 4s of particular 
interest for use in recoveries from departures and spins. A comparison of 
the stabilator power Cm. is shown in Figure 2.4. Above 60° angle-of-attack 
a full forward stick wotion (5, = +7°) has little if any effect on in- 


creasing the nose down pitching moment over 0° stabLiator setting. 


2.1.1.2 Lateral-Directional Characteristics - The lLatez:2-~ 
directional stability characteristics are improved below 40° angle of attack 
due to the addition of lesding edge slats. However, above 40° angle of 
attack the characteristics of both configurations are similar, The yawing 
moment asymmetries noted in Reference 1 on the basic F-4E are approximately 
the same as those of the siatted F-4E in both magnitude and direction, 
The directional stability characteristics at@= 0 are presented in Figure 
2.5 for 0° and -21° stabilator settings. As can be seen in Figure 2.5 
the slatted F-LE directional stability characteristics are improved over 
those of the basic F-E, 


A comparison of the lateral stability characteristics for both ean 
figurations is presented in Figure 2.6. ‘The slatted F-4E maintains a 
positive lateral stabliity throughout the angle of attack range while the 
basic F4E has almost no lateral stability in the 20° to 40° angle-of-attack 


range. 
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The addition of the leading edge elats increased both the atleron 
and the rudder effectiveness in the 20° angle-of-attack region. ‘The ad- 
verse yaw due to aileron deflection (Cys) is reduced up to a wing angle 
of attack of 30° due to the slats. These lateral-directional control 


comparisons are presented in Figures 2.7 and 2.8, 


2.1.2 Dynemic Wind Tunnel Test Results - The data obtained from 
the forced oscillation and rotary balance tests were used in generating 
dynamic derivatives for the anlytical study for the slatted F-WE, These 
dynamic derivatives are presented in Appendix A. A similar set of curves 


for the Sasic F-4E can be found in Appendix B. 


The most important of these test results is the yawing moment coef- 
ficient versus rotation rate which was obtained from the rotary balance 
tests. The leading edge slats adversely effect the yawing moment due to 
rotation rate in the 36° to 91° angle-of-attack range as shown in Figure 
2.9. In general, the propelling Cy, which is apparent on the basic F~4E 


is higher in magnitude due to the addition of leading edge slats. 


2.2 Free Flight Model Results - Free flight tests were made with a 
13% scale oe ae with and without slats. The tests were conducted in 
the NASA LEC 0 by 60 foot tunnel which included steady flight at high .. 
angles of attack and lg stalls. The same model was also dropped trom a 
helicopter. The results of these tests were used as a guide in this 
analytical study of the stall/near stall/spin regions, 


2.2.1 Free Flight Tests Results - The results from the 13% scale 
free flight model wind tunnel test are found in Reference 3. The basic F~LE 
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model was flown in a lg stall from 27° to 28" wing angle of attack before 
the model diverged in yaw. The slatted F-lE was flow to 33° wing angle 
of attack before it diverged in yav. It was noted that the use of the 
ailerons alone to control the lateral aircraft motions would cause both 


configurations to diverge at a much lower angle of attack, 


202.2 Helicopter Drop Model Results - Helicopter drops of a 13% 
scale model of the F-4E were made by NASA LRC to investigate the stall spin 
characteristics, Both the basic F-42 and the slatted F-4E were tested, 

The spin characteristics of the two aircraft are similsr. However, the 
only spin modes obtained occurred at a very low angle of attack (40° to 
50°) in sil cases. Both configurations recovered when pro-spin controls 


were neutralized, 
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3. Stall/Near Stall/Departure Characteristics 


The purpose of this section is to show the results of the stall/ 
departure analysia on the slatted F-KE aircraft and to show She changes 
in departure characteristics between the slatted and non-slatted configura~ 
tions. <A six-degree-of-freedom computer program which is described in 
Appendix C is used to calculate the stall time histories. The aerodynamic 
data used in this program is presented in Appendix A for the slatted F-HE 
and Appendix B for the non-slatted Fel4E. The production slatted F-4E weights 
and inertias are used when computing the stall/departure characteristics 
with both the sletted and non-slatted aerodynamic data. Comparisons of 
the departure characteristics due to aciodynamic changes are made for various 
laterai-directional control techniques, aircraft c.g. lecation aud stell 
entry rates. ‘Two types of maneuvers are used, an idle-power lg stall end 
a maximum power constant bank angle turn-stall. Recoveries after departure 


are also analyzed, 


Tne principle criteria used in the comparison of stall characteristics 
are angle of attack at departure and the angular acceleraiions (p und £) 
at departure. ‘The aerodynamic differences between the siatted F-hE and the 
basic F-hE cause differences ind:parture characteristics. Angle of attack 
for departure of the basic F-48 was selected as the angle of attack at which 
sideslip angle diverges since it usually departs first in yaw. Roll rate or 
yaw rate divergence is used to select the departure angle of ettack of the 


slatted F-HE since it usually departs first in roll, 


The lg stalls are started at an initial Mach of .35 and at an altitude 
of 18000 ft,, The constant bank angle turns are initiated at .6 Mech number 


at an altitude of 20000 ft.. The angle of attack is increased by inputting 
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into the computer program increasing values of airoreft+-nose-up stabilater 
as & function of time “7 tre *' - control limit (6, = <8"), A yaw guat 
was input below stall in the lg stalls to excite the roll-yaw motion, 
However the yaw gust input was not used in the constant bank angie turn 
because out-of-trim rolling and yawlng moments are already presented jn 


this type of maneuver. 


Aileron and rudder controls are used for lateral-directionally 
controlling the aircraft in the stall maneuver. The aileron is used to 
control bank angle by using a linear gain times the amount of difference 
betweer the actual bank angle and.the desired bank angle, 

5, = AlBaesirea ~ # actin) 
where A is a constant gain (A =lér 2). 
’ The rudder, when used with the aileron, controls sideslip angle as follows; 

é.. = .RB 
where Ris a constant gain (R = lor 2). 

When ‘ithe rudder control is used alone it, controls bank angle by inducing a 
email sideslip angle which creates a rolling moment due to Cy A The sideslip 
angle is also limited by the rudder to +2°, 

dy * -R, ( g ) 


desired - actual 
where R, is a constant gain (R, = 1 or 2) 


The aircraft weight and inertia configurations used in this analysis ara 
presented in Table 1. The SDP stali time histories which are presented.in 


this section of the report are summarized in Table 2. 


3.1 Effects of Flight Maneuver On Departures ~- Tho angle of attack for 
departure betxeen the slatted and basic F-4E for lg stall and constant bank 
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angle turn-stalls are summarized in Table 3 and the time histories are 
presented in Figure 3.1 through 3.6. The slatted F-4E departs at a 

higher angle of attack than the basic F-4E when maneuvering with either 
ailerons only or aileron and rudder control. The slatted F-4E does not depart 
when using the rudder only for lateral-directional control (see Figure 3.5 
and 3.6). The departure angle cf attack is lower in the constant bank angle 
turn for both the slatted and the basic F-4E at a given aileron or aileron/ 
rudder control combination. Examination of Figures 3.1 through 3.6 show the 
angular roll and yaw accelerations are much higher for the basic F=4E at 
departure than the slatted F-4E for either flight maneuver. The basic F-E 
also experiences a rapid increase in angle of attack after departure because 
of inertia coupling effects caused by the high roll and yaw rates at de- 


parture. 


3.2 Effect of Rate of Entry On Departure - The angles of attack at 


Geparture for the slatted and basic F~4E at various entry rates are yre~ 
sented in Table 4, The corresponding time histories are presented in 
Figures 3.7 through 3.14. The slatted F-lE does not depart even at the 
highest entry rate in the 1g stall maneuver while the basic F-4E depaits 

at all of the entry rates. The slatted F-4E departs at an angle of attack 
well above the basic F-4e in the constant bank angle turn (see Figures 3.12, 
3.13, and 3.14). Increasing the entry rate (@) of the basic F-4E results 
ir increasing angle of attack for departure. However increasing the entry 
rate (a) on the slatted F-lE results in ao increase in angle of atteck for 
departure since departure occurs at some finite time after the maxim 


stabilized angle of attack as determined by aircraft-nose-up stadilator. 
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As the rate of entry increases for both aircraft the amplitude of the 
laterai-directional oscillations increases, and at departure there is higher 
roll and yaw accelerstions. However, above moderately high entry rates 
(a= 5 to 10 deg/sec) the amplitude of the lateral-directional oscillations 
no longer increases or may even decrease slightly (see Figures 3.7 through 
3.14). The roll and yaw accelerations at departure also no longer increase 
with Increasing entry rate. While these trends are the same for both air- 
eraft the basic F-lE hag mach higher roll and yaw accelerations at departure 
than the slatted F-4E (see Figures 3,7 through 3.14). In addition, the 
amplitude of lateral-directional oscillations are mich higher for the basic 
F-45 then the slatted F-lE, 


3-3 Effect of Aircraft Weight, C.G., and Inertia On Departure - The 
departure angles of attack of the slatted and basic F-lE for various weight, 


@.g. and inertia configurations are presented in Table 5 (see Figures 3.1, 
3.2, 3.15, 3.16, 3.17, and 3.18 for corresponding time histories), The 
differences in departure characteristics due to aivcraft c.g. and weight 
changes are in general very minor. There isn't any consistent trend be- 
ween c.g. and the angle of attack at departure tor either the lg atall or 
the constant bank angle turn maneuvers, There is an apparent increase in 

the roll rate at departure as the c.g. was avted from forward to aft in 

the lg stalls. However, there isn’t any trend with roll rate in the constant 


bank angle turns, 


3.4 Effect of Iateral-Directional Control Inputs On Departure - The 
departure angles of attack for both aireraft using various alleron and 


rudder control techniques to maintain controlled flight are summaried in 
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Table 6 (see Figures 3.1, 3.2, 3.3 3.4, 3.5, 3.6, 3.7, and 3.12 for cor= 
responding time histories). The slatted F-4E departs at higher angles of 
attack and has lower roll and yaw accelerations at departure than the basic 
F-hp regardless of what lateral control device was used. However, the 
alleron and rudder can either improve or degrade the angle of attack at 
departure depending on the magnitude of the aileron deflection used in the 


maneuver (see Table 6). 


The basic F-4E departs at 19° (25 units) angle of attack in the lg 
stall maneuver with no lateral-directional controls while the slatted F-4E 
does not depart (see Figure 3.7). This improvement in departure char- 
acteristics is due to the increase in both lateral and directional stability 


for the slatted F-4E, 


The use of the rudder alone to control the aircraft lateral-directionally 
in the stall results in no departure for the slatted F-4E in both the lg 
stall and the constant bank angle turn-stall. The rudder is used to control 
bank angle by inducing a sideslip angle which creates a rolling moment due 
to “1g This method is very effective for the slatted F-4E because the 
dihedral effect is strong throughout the angle of attack range. The basic 
F-KE did show some improvement in angle of attack at departure over the no 


controls case in the lg stall (see Figure 3.5), 


The results of the aileron and rudder together to control the aircraft 
lateral-directionally in the stall are presented in Table 6 and show that 
departure occurs when the ailerons are used in conjunction with the rudder, 


The ailerons are limited to small deflections (6-10 degrees). The departure 
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angle of attack is high and the roll and yaw accelerations are smail for the 
Blatted F-HE, ‘The departure angle of attack (Figure 3.12) corresponds to 

the maximum attainable angle of attack using full aircraft-nose-up stabilator 
control, The basic F-4E shows a slight reduction in angle of attack at 
Geparture when the aileron and rudder controls are used in comparison to only 


rudder control. 


The ailerons used alone to control the aircraft laterally is the least 
effective of all the aileron and rudder control combinations. Departure 
angles of attack are the lowest and the angular accelerations at departure 
are the highest, While it is clear that even with the addition of slats 
the use of ailerons alone greatly restricts the angle of attack capabilities, 
however, the slatted F-4E still departs from 1° to 5° angle of attack higher 
than the basic F-4E with lower magnitudes in roll and yaw accelerations at 
departure (see Figures 3.1 and 3.2). This is due to the reduction in adverse 
yaw from the aileron and the increase in static lateral and directional 


stability of the slatted F-HE. 


3.5 Recovery From Departures - Recoveries after departure were attempted 
on both the ig stall and the constant bank angle turn maneavers. The re- 
covery method used for these calculations is full aircraft-nose-down sta- 
bilator with aileron and rudder neutralized. This ls the same procedure 
as recommend in the basic F-4E flight manual (Reference 4). ‘The time histor- 
ies showing the basic F-lE and the slatted F-4E recovery characteristics 
ere presented in Figures 3.19, 3,20, and 3.21. Both aircraft recover in 
all three cases, but the slatted F-lKE is more responsive. in the recovery 


than the basic F-hE, There is a tendency for the basic Felk to have a 
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rapid increase in angle of attack just after departure ; thereby requiring 

a longer time to reduce the angle of attack below stall compared to the 
slatted F-lHE, The reduction of roll rates and yaw rates takes a longer 
time on the basic F-4E than the slatted F-4E (see Figure 3.19 through 3.21). 
The flight manual recovery technique was quite successful on the slatted 


F-kE and it is considered to be the optimm recovery method. 


3.6 Flight Test Qualitative Results ~ In the flight test evaluation of 
the slatted F- conducted by McDonnell Douglas Corporation a large number 


of stalls were made both in lg and accelerated flight. The results of this 
test or-sram revealcd a tremendous improvement in lateral ~directional 
handling characteristics over the basic F~WE. During wings level stalls 
the angle of attack values approached 35° and above before a rolling de- 
parture would occur. In most cases departures were not experienced. At 
these high angles of attack pitch control was positive, and angle of attack 
was easily reduced by forward stick movement. Recovery was positive using 
forward stick and neutral aiilexon and rudder. Departures from accelerated 
windup turn stalls occurred in the 32° to 40° angle of attack range » and 
tended to be more gentle than the basic F-4E and are primarily in roll 
rather than yaw. Pecovery was again positive in all cases using forward 


stick and neutral aileron and rudder. 


The analytical calculations for the near stall flight conditions 
indicate the same change in departure characteristics due to the slats as 


was indicated by the flight tec’ resulta. 
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4, Spin Characteristics 
The purpose of this section is to show the results of this analysis 


and compare the basic and slatted F-4E's spin entry, spin, and spin recovery 
charcteristics at various e.g. -eroes weight configurations, The six-degree- 
of-freedom computer program which is deseribed in Appendix C is used to cal- 
culate the time histories. The serodynamic inputs used by the computer are 
presented in Appendix A for the slatted F-4E and Appendix B for the basic F-42. 
The analytical Keonetonn used for evaluating the spin characteristics is to 
start at initial lg trim conditions of .6 Mach number at 40,000 feet. A rapid 
fl. aireraft-n2se-up stabilator is input along with full right-wing-down aileron 
and fall asircraft-nose-left rudder. These control inputs are held for 39 seconds 
after which recovery controls are implemented. Recovery is accomplished with 
full aircraft nose down stabiletor, full aileron with the spin (right-wing-up) 
and rudder newbralized. The aerodynamic yawing moment asymmetries, which are 
present in the wind tunnel data above yo? angle of attack for doth the slatted 
and basic F-4E (Appendix A and B), have been applied in the same direction 
(aircraft-nose-left) as the engine gyroscopic effects so as to obtain the most 
critical. pro spin condition. Therefore all the spins evaluated in this analysis 
are to the left. The weight and c.g, configuratiors used are presented in Table 


I (see Section 3.0). A spin time history summary is presented in Table 7. 


The spin entry characteristics are compared. on the basis of the rate at 
which angle of attack and rotation rate (Q) increases with time. In the fully 
developed spin, the angle of attack and rotation rate (M) are used in comparing 
spin mode characteristics between the slatted and basic F-4E. Recovery from a 
spin is compared on the basis of the rate at which angle of attack and rota- 
tion rate (Q) decreases with time. 

4,1 Spin Entry and Spin Characteristics - The spin mode characteristics 


of the slatted F~hE is presented in Table 8. ‘The spin time histories are 
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presented in Figures 4.1 through 4.5. The effect of aerodynamic differences 
(with and without leading edge slats) on the spin entry and spin char- 
acteristics are compared at the slatted F-4E'sa weight and c.g. conditions 


(Configurations 1, 2, and 3). 


The spin times histories (with and without leading edge slats) for 
Configuration 1 (C.G, = 31.5%) are presented in Figure 4.1. The data 
without slats, has slightly higher entry rates (a and Q) than with slats 
up to time equal 16 seconés. However, after 16 seconds the rotation rate 
increases with time faster on the slatted F-4g, ‘The developed spin mode of 
both tases is considered to be a flat spin with the slatted F-4E having 
both a higher spin angle of attack as well as higher rotation rate than 


the non-slatted F-LE, 


The spin time . histories for Configuration 2 (C.G. = 29%8) ate pre- 
sented in Figure 4.2, ‘The entry rates (Gand) are about the same up 
to 20 seconds, with and without slats. After 20 seconds the slatted F-LE 
inereases in rotation rate faster with time than the non-slatted F-lE, 
The flat spin mode characteristics with and without slats are about the same 


es Configuration 1. 


For Configuration 3 (C.G. = 27%, Figure 4.3), the spin entrya@ and 
Aare sbout the same for either the slatted or non-slatted configuration. 
However after 16 seconds the slatted F-4E continues to enter ae flat spin 
(a, = 81°, Q= 130 deg/sec) but the non-slatted F-4E remains in a steep 
oscillatory spin (%, = 15° to 63°, (1= 50 to 125 deg/sec), The aerodynamic 
comparison between the slatted and non-slatted F-4E indicates that initial 


spin entry characteristics (a and{2) are about the same, but after 16 to 20 
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seconds the slatted F-4E has a greater tendency to enter the flat spin than 


the basic F-LE, 


A spin time history comparison using the aeredynamic as well as inertia 
of the slatted F-4@ (Configuration 3) and the basic F-4E (Configuration 5) is 
presented in Figure 4.4. While the inertias and gross weights of the two 
aircraft are slightly different (see Table 1) the c.g. of both configurations 
is 27¢2. The spin entry rates (& and o ) of both aircraft are about the 
Same up 16 seconds. After 16 seconds the slatted F-48 progresses much faster 
into a flat spin mode than the basic F-4E, However, unlike the previous 
case where the basic F-“E went into a .trep oscillatory spin, this time 
the basic F-lE wnt into a flat spin at the forward c.g. (279%). ‘The slatted 
F~HE has a slightly higher angle of attack and rotation rate than the basic 
F-4E in the fiat spin. The primary reason for the slightly higher angle of 
attack and rotation rates in the flet spin modes of the slatted F-HE is due 


to the higher values of propelling C, (see Section 2.1.2), 


The four spin cases that were previously discussed were deliberately 
forced into a spin by means of very rapid control inputs. It should-be 
yointed out that the initiel engle of attack entry rates { a) resulting 
from the rapid stabilator inputs were about 20 degrees per second, This 
is much higher than one would normally experience in flight. With this 
high entry rate the aircraft quickly accelerates through the angle of attack 
range (20° to Lo?) where the leading edge slats improve the lateral-direc- 
tional handling characteristics, An additional case was run (see Figure 4.5) 


to demonstrate a more normal spin entry from a constant bank angle turn using 


Configuration 1 for both aircraft. In this case no attempt is made to spin 
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the aircraft and the ailerons are used only to maintain a 65° bank angle, 
After depurture the basic F-4E progresses into a high angle of attack flat 
spin while the slatted F-~HE remains in a low angle of attack (35°-41°), 


highly oscillatory spin which is easily recoverable. 


4.2 Spin Recovery Characteristics - Anti-spin controls were applied in 
Figures 4.1 thru 4.4 at 39 seconds. ‘The recovery method of full aircraft 
nose down stabilator, aileron with, and rudder at neutral is the same as 
yecommended in the F-4E flight manual. No other recovery method was tried 
since control effectiveness of both the slatted and the basic F-4E in this 
region (a, = 70°-90°) is the same. The results show that regardless of 
the configuration if the aircraft goes into a flat spin (a, = THe to Bh°), 
recovery is mt possible with the use of existing aircraft control surfaces. 
As demonstrated by the recovery om a steep oscillatory spin (see Figure 
4.3), either aircraft recovers from any type of spin other than a high angle 


of attack flat spin. 
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5. Departure Pradiction Techniques 
Several methods have been used in the past to predict the angle cf 


attack at which an airplane becomes lateraly-directionnlly unstable thereby 
indicating the departure point of the aircraft. It is the intent of this 

section of compare several departure criteria with the results from Section 
3.0 (see Tables 3, 4, 5, and 6). ‘hese divergence criteria were calculated 


using the same aerodynamic data as used for Section 3.0, 


5.1 Directional Divergence Parameter “Dynamic Cy B - The parameter 
“dynamic Sng’ is a static derivative which includes the principle con- 


tributing terms to the undamped naturel frequency of the dutch roll mee. 
This parameter indicates directional stability with zero lateral~directional 


control deflections. ‘The expression a& described in Reference (5) is: 


Cng Dyn = Ing - (1,/%,,) “1, tana cosa>0 
The results using this parameter to predict departure are summarized in 
Figure 5.1 for both configurations, The slatted F-4B has a positive "dynamic 
Cy a through out the angle of attack range and no departures resulted during 
the SDF investigation. ‘The basic F-4s "dynamic Cn,” goes to zerc. between 
18° to 20° angle of attack. Departure ocvarred between 19° to 2° angle of 
attack during the computer investigation. This directional divergence 
parameter predicted leteral-directional co itrols-fixed departure as long as 
the entry rates tex) are low (1 to 5 degs/sec.). Generally high entry rates 
(5 deg/sec ond above) resulted in mech higher departure angles of attack 


during the computer investigation than was the predicted by Cy, B DYN values. 


5.2 Lateral Control Divergence Pardmeter - A further look at the 


lateral-directional equations of motion in terms of bank angle response 
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to aileron deflection yields another criterion. This criteria should be 

watisfied when the aileron and rudder controls are used to maintain a given 

bank angle. These expressions as diseribed in Reference (5) are ac follows: 
Aileron Alone: 


Chg, 
a—>0 Equation (1) 


“ag” “1gt 


Cnga Cnga 
Sng - Cig oa + Ky ES Cig. - ea 0) Equation (2) 
Likewise an extension of Equetion 1 for the use of proportional controls 1s 
represented by the following: 


Aileron proportional to rudder: 6, = Kod, 


Cng - cage = a J > 0 Equation (3) 
la 2 ls 

All three equations are evaluated by comparing their calculated departure 

angles of attack with those obtained from the computer time histcries of 

Section 3.0. As stated previously only the lower entry rates (@ = 1 to 5 

deg/sec) were used in these comparisons. It is felt that the departure 


angle of attack is affected very little by these lower entry rates, 


5.2.1 Aikeron Alone Divergence Parameter - The predicted angle 


of attack at departure for the aileron alone criterion (Equation 1) is 
coupared with the computer time history results (Figure 3.1, 3.2, 3.15, 
3.16, 3.17 and 3,1€) on Figure 5.2. This criterion predicted the basic 
F-4E departures within 1° to 2° angle of attack. This was not the case 
for the slatted F-lE where the criterion is unsatisfactory in predicting the 


departure angle «ttack. 
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5.2.2 Aileron/Rudder Divergence Parameter - The aileron/rudder 
criterion is calculated using Ky equal to 2. This is the same gain that 


is used for rudder proportional to sideslip in the computer results from 
Section 3.0. This divergence criterion is calculated for both the basic 
end slatted F-HE and the results are shown in Pigure 5.3. There is again 
good agreement between the SDF computer results and the departure criterion 
for the basic F-4E, however the agreement is unsatisfactory for the slatted 
F-he. ‘This criterion predicts departure well telow the actual departure 
angle of the slatted computer results. It is suspected that had K, been 

10 or larger there would have been mich better agreement with SDF results 
for the wiatted F-bE, but the basic F-4E correlation would have been un- 
satisfactory, The aileron/rudder divergence parameter is not very reli- 
eovle in particuler when strong dihedral effect is present as in the case 


of the slatted F-UE, 


9.2.3 Aileron Proportional to Rudder Divergence Criteria - 


Equation 3 is an empirical expansion of Equation 1 in Section 5.2 where 
both rudder and aileron ore consider to have a proportional relationship. 
Thus Kp is expressed as: 

by = Kab, 
The alleron and rudder are not directiy proportional to each other in the 
time histories computed in Section 3.0. The aileron is used to control 
bank angle and the rudder is used to control angle of sideslip. Thercfore 
Ko is determined by the aileron and rudder deflections at the point of 
departure (computer results, Figures 3.3, 3.4 and 3.12). These values for 
Ko are used to evaluate the alleron proportional to rudder divergence crite- 


rion (Equation 3). 
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The results of this departure criterion are summarized in Figure 


5.4, 


This criterion is unsatisfactory in predicting the angle of attack for de- 


parture particulasly for the slatted F-HE. 


5.3 f Plus 6 Axie Stability Indicator ~ A new method of determining 


tne point of lateral-directional divergence has been developed by MCAIR 


and is presented in detail in Reference 5 along with comparative free 


flight results, Because of the newness of this method a very brief deri- 


vation end explanation is appropriate, 


First we will consider the aircraft's instanteous response to sideslip 


angle about the X and Z body axes system. 


ay Srgasb 
38 I, 
a¢ _ “Igoe 
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The vector sum is called the 6 axis. As long as the A ide is in the two 
lower quadrants the aircraft is stable and no divergence is predicted. If 
the B axis should lie in the upper two quadrants or above the horizontal 

line the aircraft is unstable. The first point of instability or divergence 
ig when the 8 axis vector lies along the horizontal or along the velocity 


vector, 


In the same manner as before we will consider the instantaneous 
response to either a lateral or a directional control input or a combina- 
tion of the two. These responses are again with respect to the x and z 


body axes system, 


ad ir ed am . SLA gps 
5° “1, (8a +8) l, 
C 
ow ng, 6 4 * Chg. br qSb = Chg Sb 
6 I, (6, + &r) I, 


where 6= 89 +6r 


These vectors are shown in their positive sense in the diagram below. 


OF . Crs as 


X 
SA FS Ly 


The vector sum is called the § axis, 


When using aileron or rudder control this criteria must consider both 


the @ and the 6 axis vectors as described in the diagram on the next page. 
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This method can be applied to any control combination case or with no 


lateral-directional controls at all. 


Summaries of predicted versus actual SDF computer divergencties are 
presented on Figure 5.5 and 5.6. Four combinations of control inputs 
were looked at @long with the no control input cases. The rudder and 
aileron deflections at the angle of attack for departure (computer results) 
were used in this criterion just as in Section 5.2.3. For cases of radder 
control alone as described in Section 3.0, the rudder controls bank angle 
by producing small sideslip angles and utilizes the dihedral cffect to 
roll the sircraft, The @ and § axis method in general predicts the point 
of divergence reasonably well. ‘This method is very versatile since it is 
good for all possible combinations of aileron and rudder, aileron alone, 


rudder alone, and no lateral-diractional controls. 


MCOON 6€WELEL RPOUGLAS CORPORATION 
oe 


MDC A2181 


fp - yovEry 20 BION 


eHoad ~ 


NOTOSH MINLMVaC An-z 


(SIANSaz WEENENOD Jas) YG, 
“Ay 


200" 


cE =9 : 


SLYIS °2’T + Ghyed 


. < : ; 
(SZINSsY YZLNAKOD Jas) ~ 7 
SLYIS °S"T + ah-d LiL wos 18 
qariooo SHuALuvdsd oN ke, ae 
“ i ee se 


Be 7 x , 
/ 2ge 1/75 200 


Vat 


"Id OOCQT = “ITY SE° = WH ONES = °D°D 
STOLINDD, TYNOLIOSYLO- IVaSLYI ON HLIM XING a ks 21 ISAZY SONIM 
INS wALavesd onmvnxa 7 Ug 
ay-d Tsdon 


aayoId Wad 


S#04¢-/-O36 


Figure 5.1 


LITO IN YS A 


MDC A21U1 


Symbol Figure Maneuver 
301 1g Stall 
3.2 Turn 
3 eh 1g Stall 
3.16 lg Stall 
3.17 Turn 
3.18 Turn 
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SYMBOL 


2 where 5. = «KB 

2 : NOTE: Cpen sytbols denote basic FeLE and noldd symbols denote 
3} slatted F-4E., 
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6. Conclusions 

The following conclusions are presented based on the results of this 
investigation. 

6.1 Stali/Near Stall/Departure - The slatted F-'R shows markedly improved 
stall characteristics relative io the basic F-4E. The angle of attack at 
departure is higher in oll cases than the basic P-4E and the roll and yaw 
accelerations at devartures are mach less than with the pasic Fel. The de- 
parture characteristics of the basic F-4E usually starts with on abrupt diver- 
gence in yaw coupled with a high roll rate. The slatted F-KE usually departs 
din roll, with little yaw divergence. The evove analytically determined 
characteristics agree with the rasults obtained during flight test of the 
two configurations. © 

Application of full aircraft nose down stablintor on the slatted F-45 
departure results ina faster airplane response in reducing angle of attack. 
The basic F-LE has a repid increase in angle of attack just ufter departure 


reaulting in a longer time to reduce the angle of attack. 


Departures for the same entry rate (%) in the constant pank angle 
turn occur at a lower angle of attack than the lg stall. ‘The higher Mach 
number constant bank angle turn causes higher roll. and yaw acceleration 


for both the bagic and slatted FE. 


Increasing the entry rate (x) on beth the basic and the slattnd F-4E 
results in an increase in the angle of attack at deperture. The amplitude 
of lateral-directional oscillations increases with increasing ertry rates 
and roll and yaw accelerations are higher. However, above moderately high 
entry yates (@= 5 to 10 d@eg/sec) the oniplitude of lateral-directional 


osclilationa and the roll and yaw accelerations do not increase, 
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In summary the slatted P-45 depertures will occur at greater angles 
of attack than the basic F-4E with the rudder or rudder with restricted 
aileron contyois ond can be recovered by using Luli foward longitudinal 
control, It is also shown that with this type of alleron/rudder control 
inputs and slow entry rates that departures may not occur with the slatted 
F-lp airplane. 


6.2 Spin Characteristics - The spin mode characteristics of the basic 
P-4E and the slatted F-4E are in general the same, When nsing pro-spin 
contrels both aircraft obtained a flat spin (a, = 74°.84°) mode regardless 
of the c.g. or masa configuration. The slatted F-HE has a slightly higher 
vate of vtotation and angle of atteck in the flat spin. Recovery from the 
flat epin mode could not be made with either the basic or slatted Fur 


airplone. 


The differences in spin entry characteristius of the basic and slatted 
F-lE are smp1l when pro-spin controls are used to enter the spin, However, 
in a normal spin resulting from constant bank angle turn, the basic F-un 


is far more prone to enter a spin than the slatted F-in. 


6.3 Divergence Criteria - The departure criterion "Ch 6 ayrifunic" 
predicts the ang! > of attack for departure reasonably well. This parameter 
should only be used when the aileron and rudder sre fixed. The lateral 
control divergence parameters axe unsatisfactory in predicting departure 
engles of attack. The 8 plug & axis stability indicator shows good agree- 


ment with the couwputer time nistory departures for both intersl=directional 


controls fixed and free. 
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7, Recommendations 

(1) The Out-of-Contro2 departure recevery method of applying gixevart» 
nose-down stabilator, ailerons and rudder neutral as recanmended in the 
flight manual, is the best technique to use to recover from a departure 
with the slatted F-li aixplanc, 

(2) Recovery from a fully developed spin was studied using the computer 
reeulis, Sinee the aerodynamic characteristics of the basic F-4E ond the 
slatted F-4E are about the same in the spin regions (ap = 50° to 90°), 
the basic F-4E recommended spin recovery technique should be retained. This 
recovery method is full aircraft-nose-down stabiletor, ailerons full with 
the spin, and rudder neutralized. For a flat spin with either the basic 
or the slatted F-lE, recovery with the normal aircraft control surfaces 
4a not possible. 

(3) Sniy rudder control should be used in the high angle of attack 
region to prevent departures, However, rudder control in combination with 
restricted aileron can be used with almost the same results gas the rudder 
alone. 

(4) The anaiytical studies of a fully developed spin for the basic 
F-4E and sletted F-4E indicate little difference in the spin characteristics. 
The stall/nesr stall studies for the two aireraft indicat: the slatted F-lE 
airplane departs at higher angles of attack and is more responsive to pitch 
and lateral-directional control motions than the basic F-4E airplane. There-~ 
fore, with improved high angle of attack characteristics for the slatted F-NE 
airplanes, it 1s recommended that a stall/spin flight test program not be 


conducted on the slatted F-hs aircraft. 
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APPENDIX A 
A, Slatted F-4E Aerodynamic Data 
The data as presented in this eppendix ware used in the sis=degses- 
of-frecedom computer simulation of the slatted F-4E for both the stall/ 
Near stall and spin studies, ALL the data are input independent of Mach 
number and flexibility effects. The acrodynamic coefficients are presented, 
in the body axis system, as a function cf angle of attack ranging from . 


0° to 90°, All cata is presented at a C.G, of 33%, 


A.lL Force and Moment Wind Tunnel Data - High angle of attack static 
force and moment data on the slatted F-WE was obtained from two wind 
tunnel tests. A 5% scale model force & moment test vas conducted at NASA 
LRC in the 7 x 10 Foot High Speed Wind Tunnel. The force and moment coef- 
ficient data vere obtained at various sideslip angles from 0° to 40° angle 
of attack at .23 Mach number and a Reynolds number of 1.5 x 10° per ft. 
The 6.67% scale model force and moment test was conducted at NASA ARC in 
the 12-Foot Pressure Wind Tunnel, ‘These test data were obtained for an 
angle of attack variation from 20° to 90° at various aideslip angles for 
-23 Mach number and a Reynolds iiuibier of 4.0 x 10° per f%t.. The data from 
these two tests along with some of the data presented in Appendix B 
were used in generating the slatted F-48 aercdynawic curves presented in 


this appendix, 


The static longitudinal and lateral-directional parameters Cr, Cys 
Cms Cz, and Cy are presented in Figures A.1 through A.5, at zero degrees 


control surface deflection (5.- Sr, § 0). ‘The yawing moment coef~ 


i 


ficient C,, 1s presented in Figure A.6 as & function of sideslip angle for 


MCDONHELL DOUGLAS CORPORA?'ON 
A.3 


MCDONNELL AIRCRAFT COMPANY 2181 


APPENDIX A 


various stabilator deflections, The yawing moment asymmetries (8 = 0) of 


the wind tunnel model are presented in Figure A.7. 


The static longitudinal and lateral-directional control derivatives 
C25.) Cys.) Cag! “nga? Crea? Cos? “Ae and Ores are presented in Figure 
A.8 through A.15. The sign convention used is that a positive force or 
moment is generated by a negative value for control deflection. The control 
Gerivatives C,, and ons contain both the aileron and spoiler contribution 

0a & 


referenced to the aileron deflecticn. 


A.1.1 Adjustments in Static Pitching Moment Coefficient - Above 
fe) 
50 angle of attack the wind tunnel pitching moment coefficient curve is 
adjusted in the same manner as was done in Reference 1. That is, for a 


steady spin at a given angle of attack the following equation vias used: 


c, = teks | sinacosa (Iz - Ix) + Exz (1 - 2cos°a y] - (106) Ossina 


“Q)S¢ Equation A.1 
where the rate or rotation, angle of attack, and inertia characteristics 
were obtained from Reference 7. ‘This basic pitching moment coefficient 
is then adjusted for ¢c.g., pliching moment coefficient due to rotation, 
and stabilator deflection. The results mz; @8s = 0) are presented in 
Figure A.16. It is assumed that these pitching moment coefficient values 
obtained from the basic F-4E flight test in the spin regions would be 
applicable to the slatted F-4E, since there 1s little difference found in 


the «ind tunnel pitching moment curves between the basic and slatted F-lE, 


The pitching moment coefficients below 26° angle of attack are obtained 
by taking slatted F-HE flight test stabilator required to trim versus angle of 


attack data (Reference 8) and applying these results to wird tunnel .stabilator 


MCBONNELL DOUGLAS CORPORATION 


MCDONNELL AIRCRAFT COMPANY MDC A2181 


APPENDIX A 
pover (Cys ) to generate a pitching moment; coefficient versus angle of 
8 
attack, ‘These results along with the final pitching moment curve used 


in the study is presented in Figure A.16. 


A.1.2 Adjustments in Yawing Momen® Asyametsies - Since the ad- 


justment of yawing moment asyometries was done in the same manner for both 


the slatted and the basic F-lE, the modifications are deseribed in Appendix B. 


A.2 Forced Oscillation Wind Tunnel Data ~ Fortced oscillation testing 


was conducted in the NASA-LRCO Full Scale Tunnel on a 13% scale Model F-4E 
with leading edge slats. Three axis lateral-directional data were obtained 
from 0° to 40° angie of attack in 5° increments at a dynamié pressure of 
5.00 psf and a Reynolds number of 0,505 x 10° per f.. The frequency was 
1.0 cps and the amplitude was +5° (fy). These data along with the curves 
presented in Appendix B for the basic F-4E were used to generate the forced 
oscillation dynamic derivatives for the slatted ¥-4E, ‘The dynamic foreed 
Gscillation derivatives SL) C1» Cap? Cy? and ty. are presented in Figures. 
A.17 through A.22. The longitudinal demping derlvative Ong is the sane as 
the basic F-bE and is presented in Appendix.3. 


A.2.1 Adjustment in Roll ing Due to Roll Velocity - A comparison 
of flight test slatted Fela: (A/C 284 slatted prototype) data with SDF com- 
puter time history results shows a much higher level of roll~yaw damping 
dn the 20° to 40° angle region from the computer results. Therefore, « 
parametric study was made using a linear lateral-directional dynamie 


stability program in the 20° te 40° angle of attack region to determine 
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APPENDIX. A 
the most critical factor iavolved and the amount of change necaded to match 
flight test vesults. As expected, 1, is the most senitive parameter in 
datermining the roll-yaw damping characteristics. The magnitude of “x 
was veduced in the 20° to 40° angle of attack region to obtain more neutral 
to wistablie damping level. The original and revired curves of roll damping 


0 


are presented in Figure A.23. 


A.3 Rotery Balance Wind Tunnel Data - The rotary balanve testing was 
conticked at the NASA-LRC Vertical Spin Tunnel with a 9.1% scale Model 
F-45 with leading edge slats. Three exis componentslateralidirectional data 
were obtained from 35° to 90° angle of attack in 5° inerements at a dynamic 
pressure of 1.295 psf and a Reynolds mmber of 0.3097 x 19° per ft.. The 
axis of rotation was paralle]. to the velocity vector therefore = Oy. 
Data were obtained at epin rate coefficients (1,b/2V) of 0.0, 0.075, 
0.10, 0.15, 0.2, 0.25, and 0.30. The dynamic rotary balance coefficients 
SCn, SC], andAcy are preaentind as » function of Nyb/2v at various angles 
of attack in Figurea A.2k through A.26. ‘The longitudinal parameters 4 Cy> 


‘ 
&6,, and OC,, ave the mums as the basic F-45 and are presented in Appetty 


A.3.. Rotary Balance Data Below 35° Angle of Attack - A detail 
discussion of the method used in generating rotary balance data and ad- 


juatmenta made to the rotary balance wind tumnel data are presented in 
Appendix D. ; ; 
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B. Basic F-4E Aerodynamic Date 


The data as presented in this appendix were used in the six-degree- 
of-freedom computer simulation of the basic F-hE for both the stall and 
spin studies, All the data are input independent of Mach number and flex- 
ibility effects. The aerodynamic coefficients are in the body axils system 
as se function of angle of attack ranging from °° to go’. All data are 


presented at C.G. of 33%, 


B.1 Force and Moment Wind Tunnel Data - High angle of attack static 
force and moment data on the basic F-45 were obtained from wind tunnel 
testing, A 6.67% scale model force and moment test was conducted at NASA 
ARC in the 12-Foot Pressure Wind Tunnel, These data were cbtained with an 
angle of attack variation from 0° to go° at various sideslip angles for 
223 Mach number and a Reynolds number of 4.0 x 10° per foot. The resuits 
of this test «as used as a bases for the curves presented in this appendix 


for the basic F-4h, They are the same curves as those found in Reference 2. 


The static longitudinal and lateral-directicnal parameters Cz, Cy, Cy, 
Cy, and Cy are presented in Figures B.1 through 8.5, at zero degrees control 
surface deflection (8g, Sy; 6, = 0). The yawing moment coefficient Cy, 
is presented in Figure B.6 as a Sunotion of sideslip angle for various 
stabilator deflections, The yawing moment asymmetries (8 = 0) of the wind 


tunnel model are presented in Figure B.7. 


The static longitudinal and Lateral-directional control derivatives 
Cos 5? Cy 5» Cng.s Chea? Clea? cng. ee and Cys. are presented in 
Figures B.8 through B,15. The sign convention used is that a positive 
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force or moment is generated by a negative value for control deflection. 
The control derivatives (1s and Cs contain both the aileron and spoiler 
a 


contribution referenced to the alleron deflection, 


B.l.l Adjustment in Static Pitching Moment Coefficient - The 
basic F-LE adjusted pitching moment curve is the same as was used in Refer~ 


ence 1. 


B.le2 Adjustments in Yawing Moment Asymmetries - Both the basic 
and slatted F-lE have similar wind tunnel asymmetric yawing moments coef- 


ficient scbove 40° angle of attack. These asymmetries are positive and 
therefore promote a right spin. However, the engine gyroscopic effects 
promote a left spin. Thus the sien of the asymmetries was changed to 
create a left spin simulating a more critical condition. This is justified 
in that some aircraft will always spin to the left and some always to the 
right, which indicates that the asymmetries are aircraft dependent and 


can be applied in any manner. 


An additional modification is made to the yawing moment asymmetries 
above 70° angle of attack. This was done in order to obtain the same 
intermediate spin mode as experienced in flight test, Thus by joalenne tie 
asymmetries more positive (left spin only) the angle of attack and the 
associated rate or rotation is reduced to agree better with the results 
obtained from the basic F-48 stall/spin flight test program (Reference 7), 
The same modification is made to the slatted F-lE asymmetries. The original 
wind tunnel asymmetries along with the revised curves used in this analytical 


study are presented in Figure B,16 for both the basic and slatted F-HE, 
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B.2 Forced Oscillation Wind Tunnel Data - Yorced oscillation testing 
Was conducted in the NASA LRC 12-Foot Tunnel on a 9.1% scale Model F-4D, 
Pitch, roll, and yaw data were cbtained from 0° to 90° angle of attack 
in 5° inerements at & dynamic pressure of 5.28 psf and a Reynolds number 
of .25 x 10° per ft.. The forced oscillation Investigation was conducted 
for frequencies of 0.7 eps and 1.0 cps and angular amplitudes of 5° and 
10°, The results of this test was used as a base for the curves presented 
in this appendix for the basic F-4E, They are the same curves as those 
found in Reference 2. The dynamic forced oscillation derivatives Sing? 
i Cis Cy? Chas Cys and Cy,, are presented in Figures B.17 through 
B,23. 


B.3 Rotary Balance Wind Tunnel Data - Two separate rotary balance 
tests were conducted at NASA-LRC Vertical Spin Tunnel with a 9.1% scale 
Model F-4B. In the first test six axes data (Reference 2) were obtained 
from 35° to 90° engle of attack in 5° increments at & dynamic pressure 1.295 
psf and a Reynolds number of .3097 x 10° per ft.. The axis of rotation 
was parailel to the velocity vector, thereforeN= ee Data were obtained 
at spin coefficients (Qyb/2V) of 0.0, 0.075, 0.1, 0.25, 0.20, 0.25, and 
0,30, The axial force AC,, and pitching moment AC, rotary balance coef- 
ficients are presented in Figures B.2k and B.25. These curves are the 
same as those presented in Reference 2. The remaining coefficients have 
been changed from that presented in Refexence 2 based on the second more 
current series rotary balance tests. The updated cucfficients AC,, AC), 


AC), and ACy are presented as a function Qyb/ev at various angles of 


attack in Figures 3.26 through B.29, 
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B.3.1. Rotary Balance Data Below 35° Anle of Attack - A detail 


discussion of the method used in generating rotary balance data and edjust- 


ments made to the rotary balance wind tunnel data are presvited in Appendix D. 
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C, Description of the SDF Computer Program 
‘fhe time history date presented in this report were calculated using 
the six-degree-of-freedom computer program (SDF). The program is designed 
to determine vehicle performaiice throughout the entire flight regime of 
speed and altitude in an atmosphere and gravity field, A description of 
the basic program equations of motion and the modifications made to these 


equations for rotary balance data ere presented in the following paragraphs, 


C.l The Basic Program Equation of Motion «- The equations of motion 


as normally used in the SDF Program ere presented below. 


paved or + tee e+ Tee pq + : 
Te ae 1a ae (Equation C.1) 
eee: OE Ix, 2 ty, 2 Tenga Z 
= +2 - +x + okZ 0 | OKE Sie ATEN. a. Equation C.2 
q sana iy r eg = ee 2 (Eq ) 
I i I I 
tex - Tyra + res - ae Pr * ~fng Seng at Me (Equation ¢.3) 
z 2Z Zz a? 
: Pe : 
BH VE - Wat (Equation C.4) 
P 
vV=ewp - ur + tl (Equation ©.5) 
F F 
Ww ug- vp + = (Equation C.6) 
where . 
M, = -Tsindm sin fy Yy + Tsin Ag Cosh gy Zy (Equation C.7) 


+4 psp [ey + a Br +O), 8, + a (Cyr +o,5)] 


TeosApZy + Tein Agsin dp Xy (Equation C.8} 


sdavsé|c Paes ered ae | 
ms, 6 ay mga 


va 
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{Bquation C.9) 
rad oY (Coyx i Copp 


Cy 55 5 | (Equation €.10) 


APPENDIX ¢ 
M, = “TsinAw cosh, Xy - ToosAg Yy 
1 
+ 5 p¥sp [es + Cp, Or . Cng, 5 
’ I 2 
Fy = -Wein 9 Reon Ag + 3 pV s[¢, + 
Fy = Weos @ sin $- TsinAp cos # » + 


¥F, = Weos Scos $- Tsin Ap sin bp a 


Where Cy, Cas Spy Cx» Uys and Cy are 
and are a function of wand Bf. These 


and B. 


Parr b 

spVv Ss [&v5,8 x + oT (Cyr + C50) 
(Equatioa C.11) 

L «52 e 

zoV°S | (C, + C,, 8, + c 

. | ie a ay (eqa) | 
(Equation C.12) 

values at zero control deflections 


curves are presented in Appendix A 


C.2 Modification of Equations of Muticn For Rotary Balance Data - 
The aerodynamic characteristics of the F-4 vary nonlinearly with rotation 


rate in a spin, Since foreed oscilletion derivatives are linear it is 


necessary to use rotary balance wind 


tunnel data instead, so that the 


aerodynamic: force and moment coefficients could be obtained as a furmtion 


of rotation rate. AS described in Appendix A and B retary balance force 


and moment data were obtained by rotating the model about an axis paraliel 


to the airflow. In order to properly use these data Equations C.7 through 


C.12 must be modified. Lhe forced oscillation derivatives are applied 


only to the oscillatory components of rotation and the rotary balance data 


are applied as a Tunttion of steady state rotation, It is assumed that 


this steady state rotation is the component of rotation parallel ‘io the 


velocity vector, 
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A detaiied derivation 1s necessary to obtain the modifications to the 


equations of motion so thut rotary balance data can be used, The vector 
velationships and individual rates «ill be discussed in detail. ‘The 


total rotational velocity is given by 


N= pl +aqjy + rk (Equation €.13) 
The component off. parallel to the velocity vector is given by 
Q, =OV 
vl 
where v sul +t vj + uk 
therefore 
Q =yptygt we i 
v v {Equation C.14) 
but 
r = cos B eos a 
r = sinp 


7 = cos f sing 
Substituting into Equation C,14 gives 


Ny = peosPcosa + qsing + reosg sing (Equation C.15) 


Ae Vv 
Ay= 9 ¥ 


Dy 


The components off. along the X,¥, and 2 axes are determined as follows: 


ut 


Dy cos Boosa i +9 sing J 

+Oycos B sina k ; (Equation C.i6) 
The resulting three axes rates have been assigned the name steady state 
rotation rates and are: Pg "Dy cosB cosa 
Gg “Ny sinB 
Pg =NycosB sina (Equation C.17) 
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The rate components of the vector perpendicular to the velocity vector , 
have been assigned the name oscillatory components. The oscillatory com- 
ponents are obtained by subtracting the steady state rates Ds, Gg, and rg 


from the total rntes p, q, and re 


Po SPOS Ps 
4, = 9-4 (Equation C,18) 
Yo Fre v5 


Whom using rotary balance data, only the osc“ liutory components of 
P, Q, and + are applied to the forced oscillation derivatives. ‘The changes 
in force and moment coefficients due to steady state rotation rate are 
included as a function of M,b/2v. These rotary balance ecefficlents are 


presented as a function of a, B, ana! 


Modifying Equations C.7 through C.12 to utilize the rotary balance 
data in the manner described above gives the following results. 


My, = -fsind» sin Pp Yy + TsinAg cos Bip 2y 
+4 pVsolCy + Cy. 3p +O, Ba toe Cyx, + yp.) +AC 
lgy OF lea a ay 1y*o LpPo 1 


= TeosA 2 + Tsin A sin Pap Xyy 


d= 
H 


2am é 
* Bpv'se Ex s+ Cn * ay na) +86] 


<4 
¥ 


, = Tsindn eee fp Xy - TcosAp Yy 


2 b 
i + ice? if | 
tapv Sb |Cy + Chg pr c ba® at By on Xo + Capp,) +ACy 
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x 
F, = -Wsin @ + TeosAn + 3 pv"s [x + ox 55 5 +Ac,| 


Fy = Weos @ sin g + Tsind» cos by + Lpv's 
A 
+46, 


F, = Weos @ cos # - TainA,, sin fy + $pv*s Ic, +C,. 85 * of (Coa) ac,| 


b 
E + Cys 8x t oy Cyr, + Cy, Po) 


where 
AC , ACy, AC,» AC, ACy, and AC, are functions ofa, B, andNyb/2v and 


are presented in Appendix A and B, 
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D. Rotary Balance Data Below 35° Angle of Attack 
AS stated in Appendix A and B rotary bainnce wind tunnel data is only 


available in the 35° to 90° angle of attack range. These data ‘are necessary 
because the aerodynamic forces and moments are nonlinear with rotation rate 
in the spin angle of attack regions. Rotaky balance dats are also, needed 
below 35° angle of attack in order to ensure a smootn transition from 

below stalled flight all the way through departure. Using arbitrary 

or constant values for rotary balance data below 35° gives erratic lateral- 
directional computer results. Therefore, the method used in generating 


rotary data below 35° engle of attack is presented in the following section. 


D.L Rotary Balance Data Obtained from Forced Oscillation Derivatives - 


Rotary balance wind tunnel data is cbtained by rotating the model about 
an axis parallel to the airflow or the velocity vector. The total rotation 
rate about the velocity is called. As derived in Appendix C this rotation 
rate can be broken down into steady components as 

Ny =|Ps- a ef ss qs" 
where Ps, Ig, and dg are with respect to the x, 2, and y aircraft body axis 
system. The total rolling moment coefficient due toMy can be brokea down 
into the rolling moment cocfficients. which are dune to the individual rate 


components Dg, Xs, and qa, 3s follovs 


1 = c) +t) + Cc 
(due toNy) (due to p,} (due to rg) (due to a,) 
The individual rolling moment cemponents are the same as 
Psd 
Cc s Cy —_ 
* (ane to ps) oe 
Cc 
1 rgb 
due tor.) = ¢ 
( 3) “1, ov 
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¢ =o, ab 
lide to ay) 14 av 


Also . 
Cc Cy Dhyd 

2 2° 
(due to 2Y) Dp 2 


— 


av 


Substitution of terms yields 


Cy Nyd = ars Ped + ls psp + ‘ie qgb (Equation dD. 1) 
Nhyd ov “ov 4 “ev 
ey 


Iikewise the same method can be apolied to the total yawing moment coefficient 


and side force coefficient dueN,. 


Cy 1 psb rgd > 
ee + & + Qs 
Aye av “ Cn, ay” np “a Cn, (Equation D.2) 
2V 
C,, Nyb 
po Re, S Psd rgd asb 3 
Oyb 2 Y, “sy * ty. oy - Ney ay (Equation D.3) 
2Zv 


The relationship between the total rotation rate about the velocity vector 
Ny and its respective couwponents as derived strom Appendix C is presented 

below: 

Py, =Ny cosP cosa 

Gg =Ay *inB 

tg “Dy cosB sinc 
Assuming small values of Bthe {-JJ owing approximations are made; 

Ps “ly cosAcosa x Nycosa 

Gg =ANy sinBro {Equation D.4) 


r, Ny cos@ sinas Ny sina 
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Substj.tution of Equation D.4 into the proper Equations D.1 through D.3 


and also dividing each equation by ,b/2V yields 
Qyb 


AC = (Ci, cosa t+ Cy, ainda ) ay (Equation D.5) 
AC, = (Chy cosa + Cy sina yp . (Equation D.6) 
ACy = (Cy, cosa t+ Cy. sina ay (Equation D.7) 


where C, , Cla Cy» Cn» Cys and se are obtained from forced oscillation 
wind tunnel data. Equations 2.5 through D.7 are used to calculate constant 
slope rotary balance curves from 1° to 31° wing angle of attack. These 


curves are presented in Appendix A and B. 


D.2 Adjustment of Lateral-Directional Rotary Balance Wind Tunnel Data - 
The SDF computers time history spin characteristics indicates that the spin 
cntry characteristics of the vasic F~4E are different from the results from 
Reference 1 when the unadjusted rotury balance wind tunnel data are used 
above 35° angle of attack in conjunction with the calculated rotary balance 
data below 35° angle of attack. Since these latest results were not indicative 
of the actual airplane, the lnateral-directional rotary balance curves for the 
pasic F-LE had to be adjusted. 


As a solution to thi: problem the slopes of the yawlng aml rotiiag 
moment coefficient versus rotation rate curves are determined from both 
the generated curves (Equations D.5 and D.6) and the wind tunnel curves et 
Nyb/2v = 0. These slopes are plotted versus angle of attack and a smooth 
transition is faired between 35° and 45° (see Figure D.1). The revised 


slopes apply to the data between (]>/2V of O and +,05. This method was 
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also used to adjust the slatted F-4E rotary date (see Figure D.2). The 


final rotury belance curves are presented in Appendix A and B. A spin 
time history of the basic F-4E comparing the effect of using the updated 
rotary balance data in lieu of that used in the previous spin study 
(Reference 1) 1s presented in Figure D.3. The results show that the 
eurrent rotary balance data causes the aircraft to entry a flat spin 
slightly faster than with the data of Reference 1. This effect is attri- 
buted to the wind tunnel rotary balance yawing moment data being more 
propelling than which was used in Reference 2. The current yawing moment 
rotary balance data is compared with the corresponding previous data in 
Figure D.4 at a spin cvefficient (NQyb/2v) equal to +.05. 


D.3 The Effect of Using Rotary Balance Data Below Stall - Rotary 
balante data is used below stall. so that the same modified eauations of 


motion (see Section C.2) which utilize rotary balance data can de used in 
poth the stall/near stell/departure analysis as well as the spin analysis. 
Using rotery balance data below stall allows a continuous progrersion from 
unstalled flight through departure ali the way to the spin region. The 
purpose of this section 1s to compare the differences in departure char- 
acteristics using forced oscilistion derivatives in the normal manner (see 
Section C.1) with results obtained by applying the forced oscillation der- 
ivatives to the oscilistory components of rotation and the rotary balance 
data as a function of steady state rotation, The yawing moment derivatives 
are used as an example of the forced cscillation derivatives being applied 


in the normal manner. 


Cap By + Cp, EB = cy, (Equation D.8) 
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where p and x are the total roll rate and total yaw rate. An example of 
forced oscillation derivatives being applied to the oscillatory components 
of rotation and the rotary balance data being applied as a function of steady 


state rotation is presented telow assuming B= O, 


pee ek ER 
b 4 ro pa toe ae 
Ch, Bgh + Cy Z0P. 4 cnt) Cap 


where pp and r, are the oscillatory rate components and 


(Equation D.9) 


Po =P - Pg 

fo =r -Yg 

Ny > De” + rs” 

ee i) 4a calevlated from Equation D.6 (8 = 0) 


where Pg and r, are the stendy rate components, 


Tous combining Equations D.6 and §.9 yields 


‘b rob Q : 2 
Cp PE + Oy, TO + (Cp, cosa + Cy, sina) — = Cp,  (Equatdon 2,10) 


Vaing Equation D.4 
Ry =P “Ny cosa 
% =r -Dy sine 


Substitution of the above relations into Equation D.10 yialds 
gt La flyb rb « 0 yb 
Dp E cosa 3 Cy. E sing “ay 


+ Cap COS + Cn. csna| oe = Cue 2 Cay 


The sine and cosine terms cancel out and the result is the some as Equation D.8. 
Thus when using Equations D.5, D.6, and D.7 to generate rotary balance 


dats, the difference between using just forced osciliation data (applied to 


p and r) and forced oscillation data (applied only to py and ro} used to- 
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gether with rotary balance data (applied to p, and rg) 4s quite small. ‘the 


only approximation made was to assume that B= 0, A time history comparison 
(Figure D.5) of the two methods of applying dynamic dats shows almost 


identical trajectories up to 30° angle of attack. 


MCDONNELL DOUGLAS CORPORATION 
DS 


MOC A2161. 


= 
a 


RIC DOA RS 


aL AIRCRAFT COMPANY 


Rte & 


eee Sacre 
yee M 3 


fi 


- a 


i 
fob easken oy 
bebo Po tafe k 
t ‘ a i 


Oe mete ted 
Fi 3 pes 


ile ee ae, Sree ee 


——— 


: i: esl ' SAE Digna ey 
aa tt fb) * 6 QNW  SNOTLYNOZ ONTSN GSLVINDIVD SHOTS - 
hoo, 2 foe do Fs 2 yao TENNOL NIM SONVIVE AYVIOU HOU 
} + NoTEvIAWIS YSLAIKOD Jags NI agsn saquno 


Pan eee oe 
Hyeen dee tees “ar os ti cs ear ee 1a a0On” Dawes abes 


| 


I S80 
20 3d018- 


| ENWALOTESSOO AHO ONTMVA — 


@ONVIVE AVION JO GIOIS 


MCDONNELL eddies 1 CORPORATION 
o 


MAG 19040 IREY 29 AUG YO 


ed 


MDG A216]. 


PICDONNELL SIRCRALT COMPARY 


reas ot oer rod nese tenant rs be eemmees aad LATE 


4 
i 


SONVIVE AUVION 40 BdO1S 


ce 


“WOVE AUVION 40 Sd01S 


G ONY h SNOTLWNda ONT n aminotvs sqg01s” 

fk VIVG TSN OS GNM AONVIVE AHVIOY WOer Sa0TS" 

“ D" NOTLVINWIS “ALWHOD IGS RI ESN SHARD 0 BA07S 

sop eft 2 : ‘ vod ; a Seren rarer eer er oe i 
eee 7 YOVELY x0 mou snsuaA BEATOD waOH ve AAVOS a at 


A YSNMTOLSASOO SNEWOW: OWIAMA SL 


ee a 


D.10 


MCDONNELL DOUGLAS CORPORATICN 


MAC 10840 REY 29 AUG PO: 


piaetios 
Eo cheapean 


te | i 
ome : 


Esa ape aerenranr ent mana eens Daa 
Pp et a= |b 
oR fbeko ea aKee ae : 


} 


| aces pines Lneae oer EE 


H i B Shots 1 it 
eps fetessnrs i ara 


o ongt tie dee 4 i 


‘ ee oe ee ee a et 
cs yt 74 | 
é : i L } . 
i aan aren eee 
’ [a rig 


$e, eet aa es 
= Re eae RE OY ERR Bay 7 
Siero eae en: fees to 


pa a ee agree 


i i : a ee a a 
lgicnt i ts : ‘ot nee toga i 
Li HL ben es sn ea Se ear ; , i 

, - 7 ry 1 ¥ en 


t ee cn i i H a oor Ho . pepe] Bape 


orem os 
tigen bee 


a en ee I Sa ne inc ee 


ae as Pree gee be} 


panes eee tare ee es ice 


ea 


ria : ; shee poi jo: tt 
ae ae ar Gee rr ows i ‘6 
A ae 
Bp Nedeet op, ate nn ey PE A cary 


ae a aoa aa Fa RL 


t 


var 5 ta alerts ieee ie pee te 


bbe fob jt 


} dete se ah i beers ay, 
: : : : . st : I a i : u eet bl tod vol & : 
dae 4 


as oe et eed els eat el a jt ‘ 
west 


oe ee ee eee: ee cae LPs esac ths, i, 4a jet prs. 
. ms sy Ae. 


=a i : 


aa man ba ini 


wee ee er til tte oe Aammes  e oee st ate mee 


a oe hes ao 4a-—~ 


Me ed 


3,909 Crt ney 


Tepe ee 


of int a oo.y Re tadewede vis sik Preemie ees veep es tee 
meee zt tne : 
ie Tends nc eepater tee ele ware ied 
py - : : : . ay . ; 
OCLC ELE ELLIS L LELLILSITIS LEE SI OT MLL OPE ES SPITS mewn ng me end en wet wed 

aad : oe : , 1 
eee} er be fe Pinay ae eet ee Se wt s ee : sy ies Creek tee : 
i ae marek pane 
' Pt he ate ; ee ks, 3b begat. HS mmol 
eis i sire age ae tas Hes (aaa ; 

ia 


= 


 cicearin areata 


| ene cae Us 


- : id fe SN $i ot H ' 
Lo mare peat at ef ete ere afeceteemetecentn sfeoe foeeed see peeatiamafate foot wens 


5s eerie | 
rea 


reeeerer 
pe foe weaned me cegeeninn frmriee oak =| ea a -| art es 


(ise 
i rey 
7 a on ot a 


Fats ee (oR aTes ic OO Woes pO Hi] i ie FG Fs Se 
(ores org replete 
EE Ee pees ie te Ds Pe 
ri ese Sa a | : 
: Ga een GTR: aT | 
i 7 ‘ . 1 
t wet epee bes peered oa ho 
PS Fe et ee 
t i yee, Se : : + 
LL. a re co ote $ bef ened pet 


j rE Teen thee eet eR RE ON I 
1 ; 7 rt . ‘ x: 7 ¥ 
| ater Weer 


ey ea eae acces I $ fatrefanced abeseed oe 
aby doe 4 ‘i i | soi ‘ 
dit: ee pn 
oes oe es eer ae Hae sg ten 5 
‘aoe eae j i. iy 
poop bh caja ial ae ound areas = 
tS pe eget e OS creesepic 
ae te tieee oa ges i 
SIDER} Seaanentiens aut ieeeateded 
' oe teak I 


weet eft aa ats a Lo aa ceschta 
9 Gs foes OR AIO Gn trials Ue es fs cd oS oo 
i hot ee hy Caled 

: : : a { rye ry d 


it 


ip sense On 


. 
At. ae corse 


eure bart See Sa ta 


MOC A2181 


Fr 


RICDONNELL Ai 


CORRE, COMPANY 


i ee i 


sat eaeaael raeerres 
oe 


een ne 


\* 


ae 
“TCINGIOTSE3OO NI@S INVISNOD V 


re a tg reer oe 


2 SONMIGEN NY den YLVG hed 
— Aanis, STH Wi COVES 


Ae 7 a | Po yovaay go moe sasaan! 
A ae a eves LY THSTOLIAROD CNEOA pitts 


Figure D.4 | 
D.15 


MODONANL22. DOUGLAS CORPORATION 


MAC 10049 IREV 20 AUS YO? 


vee 


che oe] martes no = mae femmes dma feed wo Come fe no cemn het ae be momen 
Ae ae a a ee ae 
i PD ee as A Uae ae Canes | 


aia ati Ae clea ee GANT TGP 


te uk a 


in 
sepeneerttean itty o Mh 
i os) see Ae EES 


mee bee Berean 


piles Eres urea Vigas es a an send ows narra winbinice 


amen pee ee pa a a an ery nin tenner co ef aend 


ae : i : 

[ee oe aba ere eee a 
ae aot H f ; 

Seceseiieatan: vacaes dae tane aniaiar Wau aman ween! 


: i ; 
niece ‘se san eth occa aes‘ tao ‘etal: es 
tects ie ies re ; 
itt ' ry 
ee pee ee seers s. 
; he, ae : : mn 
! Be 


Prete siven s1st asa OY 


& 2 Ly 

nd ~ 
ro™ ao 3 rs 
Be ype Oe We ae 
~ hd 

& } an 
i i : i ¥ 
: 7 t , te yy ' °° 
et ete ele eee ee eit ara] Iii (eect arrae 4a, 

’ i? 5 
4 : i ue 


FATES 


omy ae mee ne neem mee ty fare _ i 


obeys 
‘ 


oo teen "2% wi gqed 7 


{ 


SC flal mentee gc rte sbesrabeas 


H * cet t ; pee tae er 
+ Tp a 
eA Re LEE LT Mi ane Lee! pees esate ts tars 


i 
ee eet ern ere 


; ee 


ar 


eee eee 
. ee 


: 4 east 
o2 cae 
- mee y 
poate Wanye 

BM Vay 


8, 


oa: 


G aoe 


no ost 


ae ne ee bee, Se en be a ee mene me 


. | ‘ aber "st, be . 724 ; 
og pee ee bg ae Pie — 


ee ns loa 


stead ation mene nm 


. 


: 

ee ee 
pit eat 
id . ° ’ Ht 


“> : ! 
Be ed 


Soe la ay Get es gs 
l CEey a 


‘ . 
wpe bee 


peaspisesicdid Sassi np isIE 
sassteny ; t nl Ae 


yee ie Shi po bebe been Paes ie i { 


e 


ta 


po anes eae aac cea oe 


tinptecee en ame vem ebie af eeeree anne 


; a “Fie ght Set ee ta ak a pry eae ees ; - : te : yee seats ra} 3 ex uay ae hagas le eres (ore Ty. 1 
“ie sjerfecied eb fed pe. pee tn ee bee vf soed ne ‘bbe ; | chad cbepenrentind e+ mfr ne 


1} 


seem io ated ot Oa Sid (ne ne ree ee ve | obo 


SER ARREST EN Se 

RS DSEAR GAT See TY 
ed 
is nett pa Stns oes, pears 


I 


ed 


t 


ep Pp 


i 
Sie off t . + t 
VOe wee Aram viene be: a 2 . a . 4 a ve . shade pmemnaes € 


pectic diwali s clescchentecnlaembenhes ta 


= e ad ae ee oean es Ns SON NOES See 
2 : ; fos 
et 
e a 3 
- Bale Diet eee bere ee 


weedy Wa tase oe abet 


: :- i ‘ 
FR EOS, a EPP TE  tS See. = 


"Or Ot 


aoe" 8 


- OF ee, 


b 


a eo Oe we 


ood 


